Fluctuation driven topological Hund insulator 
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We investigate in the framework of dynamical mean field theory a two-band Hubbard model based 
on the Bernevig-Hughes-Zhang Hamiltonian describing the quantum spin Hall (QSH) effect in HgTe 
quantum wells. In the presence of interaction, we find that a system with topologically trivial non- 
interacting parameters can be driven into a QSH phase at finite interaction strength by virtue of 
local dynamical fluctuations. For very strong interaction, the system reenters a trivial insulating 
phase by going through a Mott transition. We obtain the phase diagram of our model by direct 
calculation of the bulk topological invariant of the interacting system in terms of its single particle 
Green's function. 



Introduction - Triggered by the theoretical prediction 
[IHH and experimental observation Q of the quantum 
spin Hall (QSH) state, topological properties of non- 
interacting band-structures have been a major focus of 
condensed matter physics in recent years An ex- 

haustive enumeration of all possible topological band 
structures in the ten Altland-Zirnbauer (AZ) symmetry 
classes @ of insulators and mean field superconductors 
has been achieved using complementary approaches in 
Refs. [IMil which concludes the non-interacting clas- 
sification of these states. Since electron-electron interac- 
tions are to some extent present in every realistic system, 
the influence of electronic correlations on these topologi- 
cal bandstructures is among the key issues of this rapidly 
growing field 13 ] . 



In the presence of interactions, the Hamiltonian is no 
longer a quadratic form in the field operators and hence a 
classification based on single particle Bloch-states is not 
possible. Instead, for interacting systems, the topological 
invariant can be formulated in terms of the single particle 
Green's function (I3 - [2^ |. This was first proven 14- 20j 
for the interacting generalization of the first Chern num- 
ber characterizing the integer quantum Hall state (23. 24 1 
and has more recently been generalized to symmetry- 
protected topological states [2l], HH HH . When switch- 
ing on interactions adiabatically, the non-interacting in- 
variant connects adiabatically to the mentioned Green's 
function invariant [2lL l22|| . 



A natural question to ask is in what sense the Green's 
function topology of a system can change in the pres- 
ence of interactions. On the one hand, interactions can 
renormalize the band parameters thus giving rise to a 
new effective non-interacting band structure. On the 
other hand, dynamical quantum fluctuations can change 
the w-dcpendcnce of the Green's function which has no 
non-interacting analogue: The Matsubara Green's func- 
tion of a non-interacting system is always of the form 
G{oj,k) = (iuj-h(k))^ 1 . Deviations from this generic fre- 
quency dependence due to dynamical fluctuations have 
been shown to be capable of changing the value of the 



topological invariant by changing the pole structure of 
the inverse Green's function in w [26l423 |. Very re- 
cently, a significant simplification of the analytical form 
of the Green's function invariant has been achieved (30I — 
l32l | which makes its practical calculation much more vi- 
able. 

Main result - In this work, we study the topological 
properties of an interacting two-band Hamiltonian based 
on the Bernevig-Hughes-Zhang (BHZ) model Q for HgTe 
quantum wells. For quantum wells above a critical thick- 
ness t c , the non- interacting BHZ model shows the QSH 
effect. We consider an intra-orbital Hubbard repulsion 
U, as well as an inter-orbital term V, and a (SU(2)- 
symmetric) Hund coupling J. Then, we address the para- 
magnetic thermodynamic phase by means of dynamical 
mean field theory (DMFT). We show that starting from 
t < t c and U = V = J = 0, i.e. from a topologically 
trivial band-insulator with the low-lying "H" band com- 
pletely filled and the "E" band empty, the system can be 
driven into a QSH phase upon increasing the interaction 
strength. For the topological transition to occur it is cru- 
cial that the "Mott configuration" with equal occupation 
("H" and "E" both half-filled) is favored in the strong 
coupling regime. This happens, for instance, if a sizable 
Hund coupling J is present. The system is then driven 
away from the U = V = J — band-insulator configura- 
tion (see Fig. [2]) and the bands get closer and closer to 
one another as the interaction is increased. This effec- 
tively closes the gap and creates the conditions for the 
occurrence of topologically non-trivial phases, before the 
trivial Mott insulating state is reached at strong interac- 
tions. 

We show that, if the band gap of the non- interacting 
system is sufficiently large, the dynamical quantum fluc- 
tuations arc responsible for the topological transition, 
whereas the Hartree part of the self-energy stays in the 
trivial regime. The Mott transition at larger interaction 
strength is also associated with the low-frequency depen- 
dence of the Green's function so that its topology cannot 
be captured by static mean field theory. Our results are 
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obtained by direct calculation of the Green's function 
homotopy in the presence of inversion symmetry as pro- 
posed in Ref . [3(| . The resulting phase diagram is shown 
in Fig. H 
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FIG. 1. Phase diagram of the BHZ model with U, V = U - 
2J and J = 0.25U. For m > 2 the non-interacting model is 
topologically trivial (green) and for m < 2 is in a topologically 
nontrivial QSH phase (red). At large values of U the system 
is in a Mott phase (blue). The red empty circles at m — 
indicate that the solution is metallic. The grey lines denote 
transition regions. 

Model and methods - Previous literature on correlation 
effects in the QSH phase was mainly concerned with an 
on-site Hubbard model based on the Kane-Mele model 
for the QSH effect in graphene [33l - [36j |. This honeycomb 
model has one orbital per spin for each carbon atom, 
i.e., for each site of the honeycomb lattice. In total this 
makes two states per spin per unit cell since each unit cell 
consists of two sites. However, as an on-site Hubbard 
model, the Kane-Mele model can be seen as a spinful 
single band model. The BHZ model in contrast takes 
into account one electron-like orbital (E) and one hole- 
like orbital (H) for each spin. Hence, a lattice version 
based on this model should be naturally considered as a 
two-band Hubbard model. We consider the Hamiltonian 



H = Hbhz +Hu + H v + Hj, 



(1) 



where Hbhz is the single particle Hamiltonian as- 
sociated with the Bloch-Hamiltonian HBHz(k) = 
diag (h(k), h*(—k)), where the diagonal matrix structure 
is in spin space [37(, h{k) = (m — cos(k x ) — cos(k y )) cr z + 
Xsm(k x )a x + Xsm(ky)a y , and cr, are Pauli matri- 
ces in the E-H band pseudo spin space. The 
Coulomb interaction in a multi-orbital system can 



Hu = U^2i (nX\' n X\! + nK Q >n Xl'y an inter-orbital 
term for electrons with opposite and parallel spin, 



respectively H v = VJ2i ( n i^ n i^ + n vP n !f)j and 



Hj = (V - J)E, (n^T +n U n uJ- We also 
consider an S+iS-dike contribution, namely H$ = 

-J (ct i eti -sCp , tC nl which makes 



+" C E,t C E,i C H, J t C H,lj 

the interaction fully SU(2)-symmctric. In the follow- 
ing, we consider three forms of the interaction: i) 
H v + H v + Hj with V = U - 2J and J = Q.25U, ii) 
Hu + Hy + Hj + Hs with the same values of V and 
J and Hi) Hu alone. Hitherto, only case Hi) has been 
considered in the literature as interaction for the BHZ 
model (see Refs. 3^-41|). Yet, case in) represents a very 
unreasonable choice for a multi-orbital Coulomb vertex, 
because the inter-orbital matrix elements are entirely ne- 
glegted. Case i) and ii) are therefore much closer to a 
realistic Hamiltonian for a realizable material. 

Due to the fact that the orbitals of which Hbhz is 
constructed are not pure eigenstates of the physical spin, 
there will be a correction due to the spin-mixing influenc- 
ing cases i) and ii). However, since our results are stable 
against not too drastic variations of the parameters, we 
do expect that the spin-mixing yields only quantitative 
changes of the physics of the interacting BHZ model. 

We solve the model within DMFT (for a review see 
[42^). using the hybridization-expansion version of the 
continuous-time quantum Monte Carlo. The details of 
the code used for the numerical results can be found in 
Ref. [i^ . The advantage of this flavour of impurity solver 
is that very low temperatures can be accessed (here we 
show results for kgT= 1/100 in units of the hopping pa- 
rameter of Hbhz , but we even went down to 1 /400) and 
the calculation is sign-problem free. Since DMFT takes 
the local dynamics into account but neglects spatial fluc- 
tuations, the sclf-cncrgy depends only on the Matsubara 
frequencies: T,(iu>). The fc-dcpcndcncc of the Green's 
function comes entirely from the one-particle Hamilto- 
nian Hbhz, but the interaction changes its dynamical 
properties and hence its topology. 

We consider the half-filled case, i.e. two electrons 
on two orbitals. The E band lays 2m above the H 
band, therefore, without interaction, the system is in 
the (tih j tie ) = (2,0) configuration. Hjj alone favors 
the (1, 1) configuration. On the other hand, including 
Hv tends to restore (2,0), because inter-orbital double- 
occupation also cost energy. The Hund coupling J there- 
fore plays the crucial role for favoring (1,1) in the strong 
coupling region. The effective reduction of the energy 
difference between the orbitals is the important prereq- 
uisite to get the topological phase in the intermediate 
region. Indeed, were the bands effectively getting further 
and further apart to one another, the system would just 
forever stay in the trivial band-insulating (2, 0) configu- 
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ration. This is for instance what would happen if we only 
considered U = V and J = 0. 

In Fig. [2] we show the evolution of the orbital occupa- 
tions using cases We see how in all cases the sys- 
tem goes from the (2, 0) to the (1, 1) configuration upon 
increasing the interaction strength. Case Hi) gives an es- 
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FIG. 2. (Color online) Orbital occupation as a function of 
the interaction strength for cases i) Hu + Hv + Hj (red), ii) 
Hu + H\- + Hj + Hs + s_ (blue) and Hi) Hu (green). All plots 
show calculations performed at fcsT= 1/100, m — 3 and A = 
0.3. The vertical dashed lines mark the topological transition, 
while the Mott transition takes place when (nH,n,E) = (1, 1). 

sentially linear decrease of the orbital polarization (con- 
firming the oddness of such a choice of interaction) and, 
as a consequence, yields the topological transition in the 
weak-coupling regime where the frequency dependence 
of the self-energy is weak. As a consequence, the char- 
acter of this transition is strongly static mean-field like. 
In contrast, cases i) and ii) are much richer and physi- 
cally more interesting. Indeed, the topological transition 
is shifted to larger coupling, where the real part of the 
self-energy has a pronounced structure at low Matsubara 
frequencies. The dynamical fluctuations are therefore the 
driving force for the topological phase transition, as one 
can see in Fig. [3] where the real part of the self-energy 
for case i) is shown for different values of U. 

We now illustrate how we determine the topological 
properties from the low-frequency behavior of the real 
part of the DMFT self-energy £(wj ->• 0). The explicit 
analytical form of the topological Green's function Z2- 
invariant characterizing the QSH state as originally pro- 
posed (2H |22| employing a dimensional extension to the 
(4+l)D analogue of the quantum Hall effect [44| in the 
framework of topological field theory is rather cumber- 
some to evaluate as it involves a five-fold frequency- 
momentum integration. However, in the presence of 
inversion symmetry a tremendous simplification of this 
complicated form has recently been achieved 13011 . Fol- 
lowing also the more general analysis of Rcfs. [3ll l32j ]. 
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FIG. 3. (Color online) ^-dependence of the real part of the 
self energy. Dashed lines denote the H bands. Solid lines 
denote the E bands. All plots show calculations performed at 
inverse temperature f3 — 100, V — U — 2 J and J — U /4, m = 
3, A = 0.3. 



the recipe for the evaluation of the topological invariant 
in the presence of inversion symmetry that we employ 
here can be summarized as follows: First, one defines a 
fictitious non-interacting Hamiltonian [3~0 ] 

H(k) = -G~V = 0,k) = H(k) + S(w^0) (2) 

which has been coined topological Hamiltonian [32j ■ Ex- 
cept for values of the parameters exactly on the transition 
lines and for m = before the Mott phase, our solutions 
have insulating character and are characterized by a lin- 
early vanishing by a linearly vanishing Im S and a finite 
extrapolation of Re S(o;— >0). The latter value becomes 
large in absolute value when entering the Mott phase, as 
shown in Fig. [3] 

One can therefore proceed by calculating the Z2 invari- 
ant as proposed for the non- interac ting case in the pres- 
ence of inversion symmetry in Ref. [45| . For the reader's 
convenience, we briefly review the procedure here. First, 
one calculates the eigenstates |ri,a) of H at the four 
time reversal invariant momenta (TRIM) Ti, i = 1 ... 4, 
where a labels the occupied bands of H . Due to inversion 
symmetry, these states can also be chosen as eigenstates 
of the parity operator V with eigenvalues pi^ a = ±1. 
Next, using the phase convention y— 1 = i, one com- 
putes Si. a = ^/Pi n- Finally, the Ii invariant v can be 

5a 



expressed as 30, 



(-ir=n*. 



(3) 
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Phase diagram of the BHZ two-band Hubbard model — 
We calculate the phase diagram of Fig. [T] using the in- 
teraction ii) with U, V = U - 2 J and J = U/i. The non- 
interacting model is in the QSH phase for |m| < 2 and 
in the trivial phase for \m\ > 2. Fig. [T] shows clcary that 
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for values of 2 < m < 4, the system goes through three 
different phases as the interaction strength is increased: 
In the non-interacting limit, the system is topologically 
trivial in this regime, i.e., a normal band insulator (green 
region in Fig. [!}. At intermediate interaction strength, 
the system is driven into a QSH phase (red region) and 
in the strongly correlated limit the system is in a topo- 
logically trivial Mott phase (blue region). On the line 
separating the trivial (green) from the non-trivial (red) 
topological insulator, the system is a semi-metal. The 
Mott transition that follows is of first-order (for m = 3, 
we observe an hysteresis region of about 0.1 width in U). 
An interesting open question is therefore what happens 
for U > 12 and m > 4.5, after the two transition lines 
merge together. 

An analysis at fixed m = 3 of the evolution of the orbital 
occupation and of the self-energy shown, respectively, in 
Fig. [2] and Fig. [3j illustrates well the nature of the two 
transitions. As U is increased from the non-interacting 
case the orbital polarization gets smaller (Fig. [2]) there- 
fore the bands start to effectively get closer to one an- 
other. Since the interaction is inducing a pronounced 
frequency-structure in ReS(ioj), at [7 = 6.6 (green curves 
in Fig. [3]), the zero-frequency values exceed in absolute 
value the first eigenvalue of H , which in this case is equal 
to 1. According to Eq. ([3]) this yields a non-trivial topo- 
logical insulator. From RcS(w — > oo) we would instead 
predict a topologically trivial phase. This underlines 
the crucial role of dynamical fluctuations in the phase- 
diagram. Indeed, the high frequency limit of the self 
energy represents the static part since dynamical fluc- 
tuations decay for large absolute values of the frequency 
and only the Hartree diagrams contribute to the |w| — > oo 
limit. The Z2 = 1 phase survives till Re E(w — > 0) be- 
comes eventually larger than the largest eigenvalues of 
H (in this case equal to 5). This is the case in the Mott 
insulating state, i.e. after the occupation reaches the 
(1,1) configuration (slightly above U = 8 according to the 
red curves in Fig. [2]). Starting from a topologically triv- 
ial band-insulator, the system is therefore driven by U , 
V, and J first into a non-trivial QSH phase and finally 
into a trivial Mott insulating state. We note that the 
topological Mott insulator phase predicted in Ref. j4(| 
is qualitatively different from the Mott phase we find 
in that it describes a topological phase which occurs at 
strong interactions but can be understood at static mean 
field level where the Mott transition we observe at strong 
interactions could not take place. 

Conclusion - We studied the topological properties of 
the BHZ two-band Hubbard model with intra- and inter- 
band repulsion and a Hund coupling that makes occupa- 
tion of different orbitals with the same spin energetically 
favorable. We find that strong interaction eventually fa- 
vors the Mott state with equal orbital occupation. How- 
ever, a system with fixed trivial non-interacting band- 
parameters can go through two phase transitions the first 



of which drives the system into a nontrivial QSH phase 
and the second of which results in a trivial Mott insulat- 
ing state. The topological phase transition into the QSH 
band insulating state was shown to be driven by dynam- 
ical quantum fluctuations if the non-interacting band pa- 
rameters are chosen deep enough inside the trivial region. 
Our results stimulate the search for topological insulators 
where multi-orbital and interaction effects play a crucial 
role. This could, for instance, be the case in systems with 
d-electrons, as recently suggested for layered heterostruc- 
tures H3- 
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